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Mass  loss  from  asymptotic  giant  branch 
(AGB)  stars  is  the  most  important  driver 
for  the  evolution  of  low  to  intermediate 
mass  stars  towards  planetary  nebulae. 

It  is  also  one  of  the  most  important 
sources  of  chemical  enrichment  of  the 
interstellar  medium.  The  mass-loss  pro¬ 
cess  originates  in  the  extended  atmos¬ 
phere,  whose  structure  is  affected  by 
stellar  pulsations,  and  where  molecular 
and  dusty  layers  are  formed.  Optical 
interferometry  resolves  the  extended 
atmospheres  of  AGB  stars  and  thereby 
enables  us  to  obtain  measurements 
of  the  intensity  profile  across  this  region. 
We  present  an  overview  of  recent  results 
from  our  spectro-interferometric  obser¬ 
vations  of  AGB  stars  using  the  near- 
and  mid-infrared  instruments  AMBER 
and  MIDI  of  the  VLT  Interferometer. 


Low  to  intermediate  mass  stars,  including 
our  Sun,  evolve  to  red  giant  stars  and 
subsequently  to  asymptotic  giant  branch 
stars  after  the  hydrogen  and  helium 
supplies  in  the  core  have  been  exhausted 
by  nuclear  fusion.  An  AGB  star  is  in  the 
final  stage  of  stellar  evolution  that  is 
driven  by  nuclear  fusion,  where  a  degen¬ 
erate  carbon-oxygen  core  is  surrounded 
by  hydrogen-  and  helium-burning  layers, 
a  huge  convective  envelope  and  a  very 
extended  and  diluted  stellar  atmosphere. 
Mass  loss  becomes  increasingly  impor¬ 
tant  during  AGB  evolution,  both  for  stellar 
evolution,  and  for  the  return  of  material 
to  the  interstellar  medium.  It  reduces  the 
convective  stellar  envelope  until  the  star 
begins  to  shrink  and  evolves  at  constant 
luminosity  toward  the  hotter  post-AGB 
and  planetary  nebula  (PN)  phases,  and  is 
thus  the  most  important  driver  for  the 
further  stellar  evolution  (e.g.,  Habing  & 
Olofsson,  2003).  Mass  loss  from  AGB 
stars  is  also  one  of  the  most  important 
sources  for  the  chemical  enrichment  of 
the  interstellar  medium  and  of  galaxies. 

Depending  on  whether  or  not  carbon  has 
been  dredged  up  from  the  core  into  the 
atmosphere,  AGB  stars  appear  in  obser¬ 
vations  to  have  an  oxygen-rich  or  a  car¬ 
bon-rich  chemistry.  A  canonical  model  of 
the  mass-loss  process  has  been  devel¬ 
oped  for  the  case  of  carbon-rich  chemis¬ 
try,  where  atmospheric  carbon  dust  has 
a  sufficiently  large  opacity  to  be  radia- 
tively  accelerated  and  driven  out  of  the 
gravitational  potential  of  the  star  and 
where  it  drags  the  gas  along.  For  the 
case  of  oxygen-rich  chemistry,  the  details 
of  this  process  are  not  understood,  and 
are  currently  a  matter  of  vigorous  debate. 
Questions  remain  also  for  the  carbon-rich 
case,  such  as  regarding  the  formation 
of  the  recently  detected  oxygen-bearing 
molecule  water  in  the  inner  atmospheres 
of  carbon-rich  AGB  stars  by  the  Herschel 
space  mission  (Decin  et  al.,  2010). 

AGB  stars  experience  stellar  pulsations, 
from  semi-regular  variable  stars  (SRVs) 
on  the  early  AGB  to  large-amplitude  long- 
period  variable  stars  on  the  more  evolved 
AGB,  including  Mira  variables  and  more 
dust-enshrouded  stars  at  the  tip  of  the 
AGB.  Pulsations  and  the  induced  shock 
fronts  are  expected  to  play  a  crucial  role 
for  the  structure  of  the  extended  atmos¬ 
pheres,  but  the  details  of  these  processes 


and  their  effects  on  the  mass-loss  mech¬ 
anism  are  also  not  well  understood. 
Because  of  the  extension  of  the  atmos¬ 
phere,  temperatures  are  cool  enough 
so  that  molecules  form,  leading  to  a  sce¬ 
nario  where  molecular  layers  lie  above 
the  continuum-forming  photosphere. 

Dust  is  formed  and  believed  to  be  accel¬ 
erated  within  this  extended  atmospheric 
region  with  a  certain  condensation 
sequence,  where  dust  species  of  higher 
condensation  temperatures  form  closer 
to  the  stellar  surface  and  dust  species 
with  lower  condensation  temperatures 
form  at  larger  distances.  There  is  also  a 
notion  that  different  dust  species  may 
prevail  at  different  stages  along  the  AGB, 
along  with  increasing  stellar  luminosity 
and  increasing  mass-loss  rate. 

Optical  interferometry  at  near-IR  and  mid- 
IR  wavelengths  has  proved  to  be  a  pow¬ 
erful  tool  to  study  the  extended  atmos¬ 
phere  and  dust  condensation  zones  of 
AGB  stars,  because  of  its  ability  to  spa¬ 
tially  resolve  these  regions.  Indeed,  red 
giants  and  AGB  stars  have  historically 
been  prime  targets  for  optical  interfer¬ 
ometry,  because  of  their  brightness  and 
the  match  of  their  angular  size  to  the  typi¬ 
cal  spatial  resolution  of  optical  interfer¬ 
ometers.  However,  despite  the  long  his¬ 
tory  of  such  measurements,  they  still 
continue  to  provide  important  new  con¬ 
straints  on  the  open  questions  discussed 
above.  In  particular,  the  near-infrared 
(NIR)  instrument  AMBER  and  the  mid- 
infrared  (MIR)  instrument  MIDI  of  the  VLT 
Interferometer  (VLTI)  have  been  shown 
to  be  well  suited  for  new  measurements 
of  AGB  stars  because  of  their  unprece¬ 
dented  ability  to  provide  spectro- 
interferometric  observations  with  spectral 
resolutions  of  35-12  000  (AMBER)  and 
30-230  (MIDI).  Here,  we  report  on  very 
recent  results  (Karovicova  et  al.,  2011; 
Wittkowski  et  al.,  2011)  from  our  ongoing 
programme  to  characterise  molecular 
and  dusty  layers  of  AGB  stars  using  the 
VLTI  instruments. 


Recent  VLTI  observations 

VLTI  observations  of  AGB  stars  using 
the  instruments  MIDI  and  AMBER  started 
with  the  measurements  described  by 
Ohnaka  et  al.  (2005)  and  Wittkowski  et 
al.  (2008),  respectively.  These  studies 


24 


The  Messenger  145  -  September  201 1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

2011 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

Molecular  and  Dusty  Layers  of  Asymptotic  Giant  Branch  Stars  Studied 
with  the  VLT  Interferometer 

6.  AUTHOR(S) 


3.  DATES  COVERED 

00-00-2011  to  00-00-2011 

5a.  CONTRACT  NUMBER 
5b.  GRANT  NUMBER 
5c.  PROGRAM  ELEMENT  NUMBER 
5d.  PROIECT  NUMBER 


5e.  TASK  NUMBER 
5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

United  States  Naval  Observatory, Washington, DC, 20392  report  number 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  Messenger  145,  pgs  24-28, September  2011 

14.  ABSTRACT 

Mass  loss  from  asymptotic  giant  branch  (AGB)  stars  is  the  most  important  driver  for  the  evolution  of  low 
to  intermediate  mass  stars  towards  planetary  nebulae.  It  is  also  one  of  the  most  important  sources  of 
chemical  enrichment  of  the  interstellar  medium.  The  mass-loss  process  originates  in  the  extended 
atmosphere,  whose  structure  is  affected  by  stellar  pulsations,  and  where  molecular  and  dusty  layers  are 
formed.  Optical  interferometry  resolves  the  extended  atmospheres  of  AGB  stars  and  thereby  enables  us  to 
obtain  measurements  of  the  intensity  profile  across  this  region.  We  present  an  overview  of  recent  results 
from  our  spectro-interferometric  observations  of  AGB  stars  using  the  near-and  mid-infrared  instruments 
AMBER  and  MIDI  of  the  VLT  Interferometer. 


15.  SUBIECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

6 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


demonstrated  that  the  spectro-interfero- 
metric  capabilities  of  MIDI  and  AMBER 
are  well  suited  to  characterise  the  dusty 
and  molecular  layers  of  AGB  stars. 

We  obtained  MIDI  data  of  the  Mira  varia¬ 
ble  AGB  stars  S  Ori,  GX  Mon,  RR  Aql, 

R  Cnc,  and  X  Hya  between  2004  and 
2009.  Some  of  these  observations  were 
aimed  at  performing  a  mid-infrared  moni¬ 
toring  of  AGB  stars  in  order  to  investi¬ 
gate  intracycle  and  cycle-to-cycle  varia¬ 
bility,  and  others  were  designed  to  com¬ 
plement  AMBER  observations  to  provide 
a  more  complete  picture  of  the  dusty  and 
molecular  layers  at  the  same  time.  Mid- 
infrared  interferometric  monitoring  of  the 
sources  listed  above  has  been  performed 
by  Iva  Karovicova  during  her  PhD  thesis, 
together  with  theoretical  simulations  of 
the  expected  variability.  For  example,  the 
observations  of  RR  Aql  described  in 
Karovicova  et  al.  (2011)  include  52  obser¬ 
vations  obtained  at  13  epochs  between 
April  2004  and  July  2007  covering  three 
pulsation  cycles  (Figure  1). 

Using  the  AMBER  instrument,  we  se¬ 
cured  data  of  the  Mira  variables  R  Cnc, 

X  Hya,  W  Vel,  RW  Vel,  and  RR  Aql 
between  2008  and  2010  using  the 
medium  resolution  mode  with  a  spectral 
resolution  of  1500.  First  results  from 
these  campaigns  have  been  presented 
by  Wittkowski  et  al.  (2011). 

Some  of  these  observations  were  coor¬ 
dinated  with  concurrent  VLBA  observa¬ 
tions  of  the  SiO,  H20,  and  OH  maser 
emission  (observation  epochs  shown  in 
Figure  1),  which  provide  complementary 
information  on  the  kinematics  and  geom¬ 
etry  of  the  molecular  layers  in  which  the 
maser  emissions  originate. 


Structure  of  molecular  layers 

The  medium  resolution  (R  ~  1500)  near- 
infrared  AMBER  observations  of  the  Mira 
variables  of  our  sample  confirmed  a  char¬ 
acteristic  wavelength-dependent  shape  of 
the  visibility  function  that  is  consistent  with 
earlier  low  resolution  AMBER  observa¬ 
tions  of  S  Ori  and  that  can  be  understood 
within  the  molecular  layer  scenario.  In  this 
scenario,  the  opacity  of  molecular  layers, 
at  NIR  wavelengths,  most  importantly  H20 
and  CO,  is  larger  at  certain  wavelengths 
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Figure  1.  (Top)  Visual  light  curve  of  RR  Aql  with  our 
epochs  of  MIDI  observations  indicated  in  green. 
Epochs  of  Very  Long  Baseline  Array  (VLBA)  observa¬ 
tions  are  shown  in  blue.  From  Karovicova  et  al. 

(2011). 


Figure  2.  (Bottom)  Uniform  disc  diameter  as  a  func¬ 
tion  of  wavelength  obtained  from  VLTI/AMBER 
observations  of  the  Mira  variable  AGB  star  R  Cnc. 
The  red  line  denotes  the  best-fit  prediction  from  the 
recent  dynamical  atmosphere  model  series  CODEX. 
From  Wittkowski  et  al.  (2011). 
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so  that  the  source  appears  larger  and  the 
visibility  smaller,  and  lower  at  other  wave¬ 
lengths  so  that  the  source  appears  smaller 
and  the  visibility  larger.  Our  observations 
confirm  this  effect,  which  has  previously 
been  detected  by  interferometric  obser¬ 
vations  using  a  few  narrowband  filters. 
Hereby,  the  AMBER  instrument  has  shown 
to  be  particularly  well  suited  to  further 


study  and  to  quantification  of  this  effect 
because  of  its  wide  wavelength  coverage 
with  high  spectral  resolution.  In  our 
measurements,  the  corresponding  wave¬ 
length-dependent  uniform  disc  (UD) 
diameters  show  a  minimum  near  the 
near-continuum  bandpass  at  2.25  pm. 
They  then  increase  by  up  to  30%  toward 
the  H20  band  at  2.0  pm  and  by  up  to 
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70%  at  the  CO  bandheads  between 
2.29  pm  and  2.48  pm.  Figure  2  shows  the 
wavelength-dependent  UD  diameter  for 
the  example  of  the  source  R  Cnc. 

Recently,  new  dynamical  model  atmos¬ 
pheres  for  Mira  variables  were  developed 
that  are  based  on  self-excited  pulsation 
models  and  opacity-sampling  radiation 
treatment  (CODEX  model  series;  Ireland 
et  al.,  2008;  2011),  and  which  became 
available  at  the  time  of  our  data  analysis. 
These  models  predict  intensity  profiles 
that  are  in  excellent  agreement  with  our 
visibility  data  and  thus  the  uniform  disc 
diameters,  as  also  indicated  in  Figure  2. 

In  order  to  estimate  the  effective  temper¬ 
atures  of  the  sources  at  the  time  of  our 
observations,  we  coordinated  simultane¬ 
ous  near-infrared  photometry  at  the 
South  African  Astronomical  Observatory 
(SAAO).  The  effective  temperature  values 
obtained  are  consistent  with  those  of  the 
model  atmospheres,  as  are  the  distances 
obtained  from  the  radii  and  bolometric 
magnitudes.  Altogether,  the  agreement  of 
our  observations  with  the  CODEX  model 
atmosphere  series  increases  the  confi¬ 
dence  in  the  model  approach,  and  thus  in 
the  scenario  where  the  stratification  of 
the  extended  atmosphere  is  largely  deter¬ 
mined  by  the  pulsation  that  originates  in 
the  stellar  interior  and  in  which  the  molec¬ 
ular  layers  lie  above  the  continuum-form¬ 
ing  photosphere. 


Inhomogeneities  in  molecular  layers 

The  AMBER  closure  phase  functions  of 
our  targets  exhibit  significant  wavelength- 
dependent  non-zero  values  at  all  wave¬ 
lengths  (Figure  3  shows  the  example 
of  R  Cnc).  Non-zero  values  of  the  closure 
phase  are  indicative  of  deviations  from 
point  symmetry.  Our  data  confirm  non¬ 
zero  closure  phase  data  of  Mira  variables 
that  have  been  previously  observed  by 
other  interferometers.  Again,  the  AMBER 
data  provide  additional  information,  thanks 
to  their  spectro-interferometric  nature, 
that  helps  to  constrain  where  within  the 
atmosphere  the  asymmetric  structure 
originates.  In  fact,  we  observe  a  complex 
wavelength-dependent  closure  phase 
signal  that  correlates  with  the  features  of 
the  molecular  layers,  most  importantly 
H20  and  CO.  This  indicates  a  complex 
non-spherical  stratification  of  the  extended 


Figure  3.  Closure  phase  as  a  function  of  wavelength 
for  the  example  of  the  Mira  variable  AGB  star  R  Cnc 
obtained  with  VLTI/AMBER.  From  Wittkowski  et  al. 
(2011). 


atmosphere,  where  inhomogeneities  are 
present  at  distances  where  molecules  are 
formed  and  where  molecular  emission 
originates. 

With  the  limited  amount  of  visibility  data 
that  we  currently  have,  we  cannot  recon¬ 
struct  the  full  morphology  of  the  molecu¬ 
lar  layers.  It  would  thus  be  important  to 
conduct  an  interferometric  imaging  cam¬ 
paign  of  AGB  stars  at  high  spectral  reso¬ 
lution.  However,  since  the  visibility  data 
are  overall  well  consistent  with  spherical 
models,  we  can  already  estimate  that  the 
deviation  from  point  symmetry  originates 
from  substructure  at  a  relatively  low  flux 
level  for  an  overall  spherical  intensity  dis¬ 
tribution.  For  example,  the  values  of  R 
Cnc  in  the  water  vapour  band  shown  in 
Figure  3  could  be  explained  by  an  unre¬ 
solved  (up  to  3  milliarcseconds  [mas]) 
spot  contributing  3%  of  the  total  flux  at  a 
separation  of  about  4  mas  (for  compari¬ 
son,  the  photospheric  diameter  is  esti¬ 
mated  to  12  mas). 

There  are  several  physical  mechanisms 
that  may  lead  to  asymmetric  stellar  sur¬ 
face  structures.  For  our  targets,  we  favour 
a  scenario  where  the  pulsation  in  the  stel¬ 
lar  interior  induces  chaotic  motion  in  the 
outermost  mass  zones  of  the  extended 
atmosphere  that  leads  to  different  exten¬ 
sions  of  the  mass  zone  on  different  sides 


of  the  star,  and  thus  to  the  observed 
inhomogeneities.  An  inhomogeneous  or 
clumpy  molecular  environment  in  the 
extended  atmosphere  may  have  impor¬ 
tant  implications  for  the  non-LTE  chemis¬ 
try  in  this  region  and  may  help  to  under¬ 
stand  the  formation  of  certain  molecules 
such  as  water  in  carbon-rich  AGB  stars 
(Decin  et  al.,  2010).  It  may  also  be  the 
source  of  observed  clumpy  molecular 
features  in  the  circumstellar  environment 
at  larger  distances. 


Variability  of  dusty  layers 

The  observed  MIDI  visibility  curves  of 
RR  Aql  show  significant  wavelength 
dependence  with  a  steep  decrease  from 
8 —  9.5  pm  and  a  slow  increase  in  the 
9.5-13  pm  range.  The  shape  of  the 
8-13  pm  flux  spectrum  with  a  maximum 
near  9.8  pm  is  known  to  be  a  characteris¬ 
tic  silicate  emission  feature.  Our  rich  sam¬ 
ple  of  MIDI  data  of  RR  Aql  cover  pulsation 
phases  between  0.45  and  0.85,  i.e.  mini¬ 
mum  to  pre-maximum  phases,  for  a  total 
of  three  cycles.  For  many  different  pulsa¬ 
tion  cycles  and  phases,  we  obtained  data 
at  similar  projected  baseline  lengths  and 
position  angles.  This  gives  us  a  unique 
opportunity  to  perform  a  meaningful  direct 
comparison  of  interferometric  data  ob¬ 
tained  at  different  pulsation  phases  and 
cycles.  An  example  of  such  a  direct  visibil¬ 
ity  comparison  is  shown  in  Figure  4. 

We  concluded  that  our  data  do  not  show 
any  evidence  of  intracycle  or  cycle-to- 
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cycle  visibility  variations  within  the  probed 
range  of  pulsation  phases  (~  0.45-0.85) 
and  within  our  visibility  accuracies  of 
about  5-20%.  This  implies  either  that  the 
mid-infrared  sizes  of  the  molecular  and 
dusty  layers  of  RR  Aql  do  not  significantly 
vary  within  our  phase  coverage  or  that 
the  conducted  observations  are  not  sen¬ 
sitive  enough  to  detect  such  variations. 
The  MIDI  photometry  exhibits  a  1-2cr  sig¬ 
nature  of  intracycle  and  cycle-to-cycle 
flux  variations  which  are  most  pronounced 
toward  the  silicate  emission  feature  at 
9.8  pm.  Additional  observations  with  a 
dedicated  mid-infrared  spectrograph, 
such  as  the  instrument  VISIR  at  the  VLT, 
are  recommended  to  confirm  the  flux 
variations. 

We  performed  a  number  of  model  simu¬ 
lations  in  order  to  investigate  the  visibility 
and  photometry  variations  that  are  theo¬ 
retically  expected  in  the  8-13  pm  wave¬ 
length  range  for  the  typical  parameters  of 
RR  Aql.  We  used  a  radiative  transfer 
model  of  the  dust  shell  where  the  central 
source  is  described  by  the  dust-free 
dynamic  model  atmosphere  series,  as 
they  were  used  to  model  the  near-IR 
AMBER  data.  We  varied  model  parame¬ 
ters  such  as  the  phase  of  the  central 
atmosphere  model  and  the  opacity  and 
inner  radius  of  the  dust  shell  to  investi¬ 
gate  the  expected  photometry  and  visibil¬ 
ity  variations  during  a  pulsation  cycle. 

An  example  of  such  a  simulation  is  shown 
in  Figure  5. 

These  model  simulations  show  that 
visibility  variations  are  indeed  not 
expected  for  the  parameters  and  obser¬ 
vational  settings  of  RR  Aql  at  wave¬ 
lengths  of  8-13  pm  within  the  uncertain¬ 
ties  of  our  observations.  Variations  in 
the  flux  spectra  may  in  some  cases  just 
be  detectable.  Thus,  our  observational 
result  of  a  constant  visibility  curve  and 
only  slightly  varying  flux  spectra  at  wave¬ 
lengths  of  8-13  pm  are  consistent  with, 
and  not  contradicting,  theoretical  expec¬ 
tations  of  a  pulsating  atmosphere.  We 
used  our  simulations  as  well  to  determine 
the  best  baseline  configurations  for  the 
parameters  of  certain  targets  in  order 
to  optimise  future  interferometric  observ¬ 
ing  campaigns  that  aim  at  characterising 
the  small  expected  visibility  variations. 


Figure  4.  MIDI  visibility  amplitudes  of  RR  Aql  for 
different  pulsation  phases  (indicated  by  different 
coloured  points)  within  the  same  cycle,  chosen  to 
investigate  intracycle  visibility  variations.  From 
Karovicova  et  al.  (2011). 
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Figure  5.  An  example  of  expected  intracycle  visibility 
variation  based  on  a  simulation  using  two  atmos¬ 
phere  models  with  a  phase  difference  of  0.2  and  a 
dust  shell  that  is  located  closer  to  the  star  with  a 
higher  optical  depth  toward  stellar  minimum  com¬ 
pared  to  stellar  maximum.  The  solid  green  and  red 
lines  denote  the  two  global  models  of  this  simulation, 
the  dash-dotted  lines  the  unattenuated  stellar  con¬ 
tribution,  the  dotted  lines  the  attenuated  stellar  con¬ 
tribution  and  the  dashed  lines  the  dust  shell  contri¬ 
bution.  From  Karovicova  et  al.  (2011). 
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Characteristics  of  dusty  layers 

We  characterised  RR  Aql’s  dust  shell 
by  using  radiative  transfer  models  of 
shells  of  Al203  and/or  silicate  grains  with 
independent  inner  radii  and  opacities, 
following  earlier  such  attempts  in  the  liter¬ 
ature.  We  obtained  best-fit  results  for 
RR  Aql  with  a  silicate  shell  alone.  The 
addition  of  an  Al203  shell  did  not  result  in 
any  improvement  in  the  model  fits.  The 
best-fit  model  for  our  average  pulsation 
phase  of  &v  =  0.64  ±  0.15  includes  a  sili¬ 
cate  dust  shell  with  an  optical  depth  of 
tv  =  2.8  ±  0.8  and  an  inner  radius  of 
Rin  =  4.1  ±  0.7  RPhot  (where  RPhot  is  the 
photospheric  radius)  and  uses  a  central 
intensity  profile  corresponding  to  an 
atmosphere  model  with  an  effective  tem¬ 
perature  of  2550  K.  Figure  6  shows  a 
comparison  of  our  radiative  transfer  model 
to  the  RR  Aql  visibility  data  of  one  of 
our  epochs.  We  conclude  that  a  radiative 
transfer  model  of  the  circumstellar  dust 
shell  that  uses  dynamical  model  atmos¬ 
pheres  representing  the  central  stellar 
source  can  reproduce  well  the  spectral 
shape  of  both  the  visibility  and  the  pho¬ 
tometry  data  (not  shown). 

In  addition  to  RR  Aql,  we  determined 
best-fit  dust  shell  parameters  for  other 
sources  of  our  sample  (S  Ori,  GX  Mon, 
and  R  Cnc).  The  mid-IR  visibility  and 
spectra  of  each  star  can  be  described  by 
an  Al203  dust  shell  with  a  typical  inner 
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Figure  6.  Radiative  transfer  model  of  RR  Aql  (solid 
lines)  compared  to  the  visibility  data  of  one  of  our 
epochs  including  two  observations.  The  solid  lines 
indicate  our  best-fit  model,  while  the  contributions  of 
the  attenuated  stellar  and  dust  components  alone 
are  denoted  by  the  dotted  and  dashed  lines,  respec¬ 
tively.  From  Karovicova  et  al.  (2011). 

radius  of  2-2.5  stellar  radii  and/or  a  sili¬ 
cate  dust  shell  with  a  typical  inner  radius 
of  4-5  stellar  radii.  This  result  is  con¬ 
sistent  with  the  scenario  of  a  dust  con¬ 
densation  sequence  where  Al203  grains 
form  closer  to  the  stellar  surface  and 
silicate  grains  at  larger  radii,  which  is  as 
well  consistent  with  the  condensation 
temperatures.  Our  best-fit  values  of  the 
optical  depths  of  the  Al203  and  silicate 
dust  shells  together  with  the  mass  loss 
rates  adopted  from  the  literature  may 


confirm  earlier  suggestions  of  a  sequence 
where  the  dust  content  of  stars  with  low 
mass-loss  rates  is  dominated  by  Al203 
grains,  while  the  dust  content  of  stars 
with  high  mass-loss  rates  predominantly 
exhibit  substantial  amounts  of  silicates. 
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A  composite  of  images  from  VISIR  and  NACO  of  the 
red  supergiant  star  Betelgeuse.  The  image  is  5.63 
by  5.63  arcseconds  in  extent  and  shows  the  exten¬ 
sion  of  the  nebula  around  the  star  to  be  some  26 
milli-parsec  (550  AU).  The  supergiant  star  itself  is 
resolved  by  interferometry  at  44  milli-arcseconds 
(3.4  AU).  The  larger  image  is  taken  with  VISIR  in  ther¬ 
mal-infrared  filters  sensitive  to  dust  emission  and 
the  compact  image  (inside  the  occulting  spot)  is 
taken  with  NACO  in  J-,  H-  and  K-band  filters.  Some 
of  the  non-spherically  symmetric  eruptions  close  to 
the  star  in  the  NACO  image  can  be  traced  to  much 
larger  distances  in  the  VISIR  image.  More  details  can 
be  found  in  Kervella  et  al.,  A&A,  531,  117,  2011  and 
in  Release  eso1121. 
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